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Renormalizations,  which  sum  classes  of  diagrams  to  all  orders  in 
perturbation  theory,  are  applied  to  the  electron  propagator.  The  sepa- 
ration of  relaxation  and  correlation  self-energy  diagrams  enables  the 
construction  of  electron  propagators  that  include  relaxation  contribu- 
tions to  all  orders.  Correlation  self-energies  are  obtained  in  second 
order,  third  order,  and  in  the  diagonal  2p-h  Tamm-Dancoff  approximation, 
which  is  also  a  renormalization.  These  propagators  are  applied  to  the 
electron  binding  energies  of  CN~,  NFL",  PFL~  and  Cl2-  An  alternative 
propagator  renormalization  through  the  transition  operator  method 
reference  state  is  applied  to  the  electron  binding  energies  of  Ne,  FLO, 
NH2"  and  CN". 

The  anti symmetrized  geminal  power  wavefunction,  which  is  the  form 
of  the  consistent  ground  state  of  the  random  phase  approximation  of  the 
polarization  propagator,  is  variational ly  optimized  for  the  LiH  ground 


vm 


state  potential  curve.  Random  phase  approximation  polarization 
propagator  calculations  yield  excitation  energies. 

Physical  concepts  suggested  by  these  theories  are  discussed, 


IX 


INTRODUCTION 

Modern  chemical  investigation  inevitably  revolves  about  the  axis  of 
molecular  electronic  structure  theory.  The  quantum  mechanical  investiga- 
tion of  the  structure  and  properties  of  isolated  molecules  lays  a  founda- 
tion for  studies  of  complex  chemical  systems.  Energy  differences  between 
states  are  properties  of  continuing  importance  in  molecular  electronic 
structure  theory,  for  they  provide  insight  into  macroscopic  chemical 
phenomena.  Propagator  theory  is  a  theoretical  technique  that  concen- 
trates on  determining  energy  differences  and  transition  probabilities. 
The  concepts  that  are  introduced  with  propagator  theory  therefore  explain 
the  energy  differences  themselves  without  necessarily  treating  the  total 
energies  or  wavefunctions  of  the  states  of  interest. 

Molecular  energy  differences  are  associated  with  nuclear  and  elec- 
tronic rearrangements.  The  Born-Oppenheimer  approximation  permits  a 
separation  of  electronic  and  nuclear  motion  that  facilitates  quantum 
mechanical  calculations  based  on  the  time  independent  Schrodinger 
equation.  With  the  nuclei  fixed  at  given  positions,  the  electronic 
states  of  the  molecule  are  studied.  A  series  of  nuclear  configurations 
defines  an  adiabatic  potential  energy  surface,  where  the  total  electronic 
energy  plus  the  internuclear  repulsion  energy  is  plotted  against  the 
nuclear  positions.  Ground  state  potential  energy  surfaces  have  typically 
been  calculated  using  the  variational  method;  propagator  theory  has  made 
few  contributions  to  this  field.  However,  for  the  calculation  of  energy 
differences  between  potential  energy  surfaces,  propagator  theory  is 
competitive  with  current  variational  techniques.  In  the  chapters  that 
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follow,  propagator  theory  will  be  applied  to  energy  differences  between 
electronic  states.  The  calculations  will  be  made  in  a  limited  basis  of 
square  integrable  functions  and  there  will  be  no  explicit  treatment  of 
the  continuum. 

Electron  binding  energies  are  energy  differences  between  N  electron 
and  N±l  electron  states,  i.e.,  ionization  energies  and  electron  affini- 
ties. Hartree-Fock  theory  provides  predictions  of  electron  binding 
energies  based  on  Koopmans's  theorem  (Koopmans,  1933)  which  states  that 
orbital  energies  are  interpretable  as  ionization  energies  or  electron 
affinities.  Unfortunately,  the  invariance  of  the  Hartree-Fock  wavefunc- 
tion  under  unitary  transformations  of  the  occupied  (or  the  unoccupied) 
orbitals  among  themselves  places  in  doubt  the  validity  of  the  orbital 
energy  interpretation.  Propagator  theory  justifies  and  delimits  the 
simple  orbital  model  and  suggests  ways  to  improve  it.  These  ways  include 
the  extension  of  the  orbital  concept  to  quasi-particle  concepts  of  elec- 
tron binding  energies  and  the  introduction  of  energy-dependent  correla- 
tion potentials. 

Excitation  energies,  or  energy  differences  between  electronic  states 
with  the  same  number  of  electrons,  are  also  amenable  to  propagator  inter- 
pretations. Hartree-Fock  theory  also  gives  a  model  for  excitation 
energies  that  can  be  extended  with  perturbative  propagator  corrections, 
but  propagator  theory  also  suggests  models  that  abandon  the  orbital 
picture  as  a  point  of  reference.  Correlated  model  Hamiltonians  provide 
physical  pictures  for  the  ground  and  excited  states  and  explain  the 
nature  of  the  interactions  present  in  each  state. 

These  applications  illustrate  how  propagators  provide  direct  deter- 
mination of  energy  differences  and  introduce  concepts  that  are,  so  far, 
inaccessible  to  purely  variational  theories. 


CHAPTER  ONE 
ELECTRON  BINDING  ENERGIES 

1.1  The  electron  propagator 

Electron  propagators  efficiently  determine  and  explain  electron 
binding  energies.  The  poles  of  the  electron  propagator  are  energy 
differences  between  the  N  electron  ground  state  and  the  N±l  electron 
states,  i.e.,  electron  binding  energies.  The  residues  of  the  propagator 
are  overlap  amplitudes  between  the  N  electron  and  N±l  electron  states 
and  describe  the  electronic  structure  of  the  ground  state.  The  electron 
propagator  matrix  (Linderberg  and  Ohrn,  1973),  G,  in  a  spin-orbital 
basis  has  the  elements 


fn*(i)fn(j)       gn*(i)gn(J) 

GiJ(E)  =  I   E-En(N+l)+E0(N)+ia+  \   E+En(N-l)-EQ(N)-ia  '       (1_1) 


The  total  energy  of  the  N  electron  ground  state  is  EQ(N),  and  E  (N±l)  is 
the  total  energy  of  the  n-th  N±l  electron  state.  The  poles  of  the 
propagator  are  real  energy  differences  when  the  positive  parameter  a 
approaches  zero.  The  continuum  can  be  treated  with  the  electron  propaga- 
tor in  conjunction  with  other  techniques  (Mishra,  1981).  The  overlap 
amplitudes  (Purvis  and  Ohrn,  1975)  between  the  ground  and  N±l  electron 
states  are 

fn(i)  =  <N+l,n|ai+|N>  gn(1)  =  <N-l,n|ai|N>  ,    (1-2) 

where  the  field  operators  satisfy  the  relations 


[a,,aj]+  .  [a^.a/j,  .  [.,,«  ^  -  «   .  (1-3; 


In  this  notation,  the  many-electron  Hamiltonian  is 


H    -     H   h.-aJa.  +  k  I  H  I   <ij||kl>a.+a.+a  a  ,        (1-4; 
i  j  1J  1  J    i  j  k  1        i  J  i  k 


where  <ij||kl>  =  <i(l)j(2) 


1-P12 


rr2' 


k(l)l(2)> 


NUCLEI   Z 

hf1  ■  <i|-^v2+  I       — -a-|j>  .  (1-5; 

J  a     R  -r 


This  Hamiltonian  operates  in  Fock  space  and  is  independent  of  the  number 
of  electrons  in  the  states  on  which  it  operates. 

1.2  An  uncorrected  propagator 

Consider  the  model  Hamiltonian 

H0  -  I  eiai+ai  (1-6) 

in  the  spin  orbital  basis  that  diagonal izes  the  Fock  matrix, 

Fij     -     hij  +  n<iMm><akV   .  (1-7) 


where 


t  t 

<ak  a.j>  =  <ground  state |ak  a-,  | ground  state>  ,  (1-8] 


so  that 


Fij  ■  «i«u  •  (i-s: 


Since  the  model  Hamiltonian  contains  only  one-electron  terms,  the  elec- 
trons do  not  interact.  The  model  ground  state  is 

N    t,  , 

|N>  =   n  a.'|vac>  (1-10) 

i=l  n 


where 

N  |vac>  =  I   n. |vac>  =  I   a,  a-|vac>  =  O|vao  .       (1-11) 


The  ground  state  total  energy  is 


<N|HQ|N>  =  I   e.O^a^  =  J  ^n.  .  (1-12) 


There  are  N  orbital s  where  n.  =  1;  the  rest  of  the  orbital s  have  n-  =0. 
The  N±l  electron  states  differ  from  the  ground  state  in  the  choice  of 
occupation  numbers.  The  poles  of  the  model  electron  propagator  are 

EP0LE  =  EoW  "  M^  =  en  (1"13) 

which  are  negative  ionization  energies,  and 

EP0LE  =  En^N+1)  "  E0^  =  en  (1"14) 

which  are  negative  electron  affinities.  The  overlap  amplitudes  are 

f  (i)  =  <N+l,n|a.+|N>  =  6. 
n  1        in  (1-15) 

gn(1)  =  <N-l,n|a.|N>  »  &-n   . 

The  purpose  of  this  exercise  is  to  recover  Koopmans's  theorem 
(Koopmans,  1933)  and  the  attendant  orbital  concept  of  electronic 


structure.  The  remaining  discussion  considers  how  the  electron  propagator 
can  be  extended  to  describe  correlated  electron  dynamics. 

1 .3  Toward  correlated  electron  propagators 
Given  the  notation  for  the  propagator 


G.j(E)  =  «ai;aj-+»E  ,  (1-16; 


the  identities 


,  E   (N+l)-En(N) 

E[E  -  En(N+l)   +  EqCN)]1     -     1  +  E-En(N+l)+EQ(N)    • 


(1-17) 


<N|a.|N+l,n>(En(N+l)   -  EQ(N))     -     <N|[a.,H]|N+l,n>   , 
and  the  corresponding  relations  for  the  N-l  electron  case  imply  that 


E«ai,ajt»£     =     <[ai,aj+]+>  +  «[ai,H];aj+»E   .  (1-18) 


This  equation  can  be  iterated  to  give 

«ai,aJ.+»E  =  E'^Ca-.a^1"]^  +  E"2<[[ai  ,H],aj+]+>         (1-19) 
+  E"3<[[[a.,H],H],aJ.+]+>  +  ...  . 

The  definition  of  superoperators  (Goscinski  and  Lukman,  1970),  t  and  H 

IX  =  X  HX  -  [X,H]  ,  (1-20) 

which  act  on  the  set  of  field  operators  (Redmon  et  al . ,  1975) 


X  =  {a^a.+a^a.+a/a^,...}  ,  (i-2i; 


and  the  definition  of  a  scalar  product  (Goscinski  and  Lukman,  1970) 

(X|Y)  =  <[Y,X+]+>  (1-22) 

convert  the  expansion  in  equation  1-19  to  the  form 

«ai;a-+»E  =  (a^EI-H)"^)  .  (1-23) 

The  inverse  operator  which  occurs  in  the  propagator  matrix, 

G(E)  =  (a|(EI-H)_1a)  ,  (1-24) 

can  be  circumvented  by  the  inner  projection  technique  (Lowdin,  1965): 

G(E)  =  (a|h)(h|(EI-H)h)_1(h|a)  ,  (1-25) 

where  the  operator  manifold  is  complete,  i.e., 

h  =  {a, a  aa^a  aaa,...}  ,  (1-26) 

and  where  the  superoperator  matrix  elements  are  evaluated  with  respect  to 
the  exact  ground  state.  A  partitioning  (Lowdin,  1965  and  1971)  of  h, 

h  =  {a|f}  ,  (1-27) 

such  that 

(a|f)  =  0  (a|a)  =  Cf if)  =  1  »  (1-28) 

implies  that 


G_1(E)  =  (a|(Et-H)a)  -  (a|Hf )(f | (EI-H)f )_1(f |Ha)  .        (1-29) 


For  the  model  uncorrelated  system  where 


|N>  =  |Hartree-Fock>  (1-30) 

and  where  the  model  Hamil tonian  is 


H0  =  I   eiaiTai  (1-31] 


the  exact  inverse  propagator  is 

G0_1(E)  =  (a|(EI-HQ)a)  ,  (1-32) 

since  the  second  term  in  the  right-hand  side  of  equation  1-29  vanishes. 
Koopmans's  theorem  appears  once  again  since 

[^^(E)]^  =  (ai|(Et-H0)aj)  =  (E-e^j  .  (1-33) 

For  the  full  Hamil tonian,  approximations  of  the  exact  ground  state  and 
choices  of  a  subspace  of  f  are  necessary.  If  the  ground  state  is  written 
as 

|N>  =  |Hartree-Fock>  +  |correlation>  ,  (1-34) 

then  the  propagator  becomes 

G_1(E)  =  El  -  (a|Ha)Hp  -  (a|Ha)CQRR  (1-35) 

-  (a|Hf)(f|(EI-H)f)_1(f|Ha)  , 

« 

where  (a|Ha)HF  is  (a|Ha)  taken  only  with  respect  to  the  Hartree-Fock 
contribution  to  the  ground  state,  and  where  the  remainder  of  (a|Ha)  is 
(a|Ha)£QRR.  In  other  words, 

G-X(E)  =  Gq-^E)  -  I'    -  l'(E)    .  (1-36) 
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The  superoperator  matrix  element 

(ajfia.)     =     h..   +  I  I  ^kMilxa+a^  (1-37) 

depends  only  on  the  one-electron  reduced  density  matrix  and  is  indepen- 
dent of  E.  Taken  together,  the  energy  dependent  and  energy  independent 
corrections  to  the  uncorrected  propagator, 

V   +X'(E)  -     1(E)  (1-38) 

are  the  self -energy  and  appear  in  the  Dyson  equation, 

G_1(E)  =  G0_1(E)  -  1(E)    .  (1-39) 

The  self-energy  contains  all   correlation  effects. 

1.4  Approximations 

Choices  for  the  ground  state  and  for  a  subset  of  f  define  approxi- 
mate propagators.  A  specific  physical  problem  or  computational  feasibil- 
ity may  motivate  these  choices.  Redmon  et  al .  (1975)  and  Nehrkorn  et 
al .  (1976)  thoroughly  analyzed  the  corrections  to  Hartree-Fock  from 
single  and  double  excitation  contributions  to  the  second  order  Rayleigh- 
Schrodinger  perturbation  theory  wavefunction.  Several  authors  have 
studied  choices  for  the  operator  manifold  (Pickup  and  Goscinski,  1973; 
Purvis  and  Ohrn,  1975;  Redmon  et  al.,  1975;  Nehrkorn  et  al . ,  1976; 
Herman  et  al . ,  1980a  and  1980b;  Simons,  1977  and  1978;  McCurdy  et  al . , 
1977).  An  example  is  f.,: 

U     =     {fklm}  =  {(ak+alam)  +  <nk>al6km  "  <nk>am6lm)Nklm^}  ' 
Nklm  =  <nk>(1  "  <V  "  <nm>}  +  <nl><nm>  •  ^~^ 
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The  components  of  the  self-energy  also  derive  from  diagrammatic 
perturbation  theory  (Cederbaum  and  Domcke,  1977).  In  this  approach,  the 
terms  of  the  self-energy  are  classified  by  their  order  of  electron 
interaction.  Born  has  discussed  the  connections  between  the  order  by 
order  approach  and  the  inner  projection  formalism  (Born,  1979).  Operator 
product  manifolds,  such  as  f.,  or  fr,  sum  classes  of  terms  in  the  diagram- 
matic expansion  to  all  orders  in  electron  interaction.  The  infinite 
summations,  or  renormalizations,  are  useful  when  the  convergence  of  the 
diagrammatic  series  is  in  doubt.  For  example,  the  choice  of  the 
Hartree-Fock  ground  state  and  the  f,  operator  product  manifold  sums  all 
terms  in  second  order  and  also  includes  classes  of  diagrams  called  rings, 
ladders,  and  mixed  ring-ladder  diagrams  to  all  orders.  If  the  ground 
state  includes  double  excitations  from  the  first  order  RSPT  wavefunction 
and  single  excitations  from  the  second  order  RSPT  wavefunction,  then  the 
renormalization  sums  all  diagrams  through  third  order,  in  addition  to 
other  diagrams  in  higher  orders.  The  two  complementary  modes  of  analysis 
of  the  self-energy,  from  diagrammatic  theory  or  from  the  approximation 
of  the  inner  projection  formula,  enable  verification  of  approximate  self- 
energies.  They  are  alternative  terminologies  which  suggest  improvements 
of  the  self-energy. 

Yet  another  method  of  analysis  of  the  self-energy  is  the  identifica- 
tion of  so-called  relaxation  and  correlation  terms  by  perturbation  theory 
(Born  et  al.,  1978).  Relaxation  terms  are  terms  that  occur  in  the 
difference  between  Hartree-Fock  total  energies  for  the  N  electron  and 
N±l  electron  states.  (The  ASCF  method  is  discussed  in  Bagus,  1965.)  The 
perturbation  is  the  difference  between  the  Fock  operators  of  each  state. 
In  general,  the  total  energies  are  unrestricted  Hartree-Fock  total 
energies.  The  relaxation  terms  appear  in  the  diagonal  elements  of  the 
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propagator  matrix  when  E  =  e . ,  where  e-   is  the  orbital  energy  pertaining 
to  the  N±l  electron  state  of  interest.  Perturbation  theory  connects  the 
SCF  total  energy  of  the  N±l  electron  states  with  the  orbitals  and  orbital 
energies  of  the  N  electron  SCF  calculation.  The  terms  in  the  self-energy 
which  are  not  relaxation  terms  are  correlation  terms.  Since  all  self- 
energy  contributions  treat  correlation  effects  with  respect  to  the 
uncorrected  Hartree-Fock  propagator,  the  use  of  the  term  correlation 
self-energy  will  refer  only  to  contributions  that  do  not  originate  from 
relaxation  effects. 

A  further  approximation  proceeds  from  investigations  (Schirmer  et 
al.,  1978)  which  demonstrate  the  unimportance  of  nondiagonal  elements  of 
the  propagator  matrix  when  the  matrix  is  constructed  in  the  basis  of 
canonical  Hartree-Fock  orbitals.  (This  observation  holds  for  outer 
valence  and  core  binding  energies  of  molecules  near  their  equilibrium 
geometries.)  Neglecting  the  nondiagonal  elements  sacrifices  little 
numerical  accuracy.  Mixings  of  canonical  Hartree-Fock  orbitals  of  the 
same  symmetry  in  the  overlap  amplitudes  are  obliterated  when  the  diagonal, 
or  quasi-particle,  approximation  is  made.  Previous  studies  of  the  non- 
diagonal  propagator  have  also  indicated  that  these  mixings  are  small 
(Schirmer  et  al.,  1978).  The  quasi-particle  approximation  causes  many 
poles  to  be  assigned  to  each  canonical  orbital.  In  the  uncorrected 
propagator,  only  one  pole  is  assigned  to  each  orbital. 

1.5  Approximate  propagator  expressions 

With  the  orbital  energy  and  the  relaxation  terms  in  the  self -energy 
included  in  the  ASCF  energy  difference,  the  renormalized  quasi-particle 
propagators  under  consideration  will  have  the  form  (Kurtz  and  Ohrn,  1978) 
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G'^U)  =  E  -  (ASCF)i  -  IC0RR..(E)  .  (l-4i; 


The  quasi-particle  assumption  permits  direct  iteration  of  equation  1-41 

with  respect  to  E,  i.e.,  there  is  no  need  to  diagonal ize  with  every   new 

value  of  E.  However,  since  relaxation  terms  are  strictly  defined  only 

when  E  =  e • ,  the  expression  (Kurtz  and  Ohrn,  1978) 

G'^U-j)  =  e1  -  (ASCF)i  -  ^^-(e.)  (1-42) 

will  be  evaluated  and  reported.  Typically,  when  the  quasi-particle 
approximation  is  valid,  dG~  (E)/dE  is  small.  Consequently,  the  results 
of  1-41  and  1-42  closely  agree  (within  0.1  electon  volts).  Diagrammatic 
theory  suggests  the  second  and  third  order  correlation  self-energies. 
(See  the  Appendix  for  correlation  self-energy  formulas.)  Another  approx- 
imation, which  derives  from  the  inner  projection  formalism  with  the  f- 
operator  product  manifold  and  a  Hartree-Fock  ground  state,  ignores  off 
diagonal  terms  in  the  matrix  (Born  and  Ohrn,  1978,  1979  and  1980) 

(f3|(EI-H)f3)"1  .  (1-43) 

This  approximation  is  the  so-called  diagonal  2p-h  Tamm-Dancoff  approxima- 
tion (2p-h  TDA)  (Cederbaum  and  Domcke,  1977;  Schirmer  and  Cederbaum,  1978) 
and  is  no  more  difficult  to  compute  than  second  order.  It  renormalizes 
certain  ring,  ladder  and  mixed  ring-ladder  diagrams.  (For  a  discussion 
of  the  nondiagonal  2p-h  TDA,  see  Mishra  and  Ohrn,  1980.)  All  of  the 
following  applications  of  the  ASCF  plus  correlation  self-energy  propaga- 
tors are  discussed  in  Ortiz  and  Ohrn,  1980  and  1981. 
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1.6  Vertical  and  adiabatic  electron  affinities  of  CU 

A  vertical  electron  binding  energy  is  an  energy  difference  between  N 
electron  and  N±l  electron  states  at  the  same  nuclear  geometry.  An 
adiabatic  electron  binding  energy  is  an  energy  difference  between  the 
same  states  evaluated  at  their  equilibrium  nuclear  geometries.  For  Cl2, 
the  change  in  the  eauilibrium  internuclear  distance  upon  the  addition  of 
an  electron  is  about  1.3  bohr.  To  calculate  the  adiabatic  electron 
affinity  of  Cl?  requires  the  difference  between  the  energies  at  the 
internuclear  distances  of  the  neutral  molecule  and  of  the  anion.  Either 
Cl?  or  Cl2"  curves,  combined  with  the  appropriate  vertical  electron 
affinity,  give  the  adiabatic  electron  affinity. 

Basis  sets  for  all  the  systems  consist  of  contracted  Gaussian  type 
orbitals  (CGTOs).   Dunning  and  Hay's  (1977)  chlorine  basis  set  is  the 
simplest  choice  (Basis  A  in  Table  1-1);  d  functions  and  diffuse  s  and  p 
functions  are  subsequent  augmentations.  As  the  basis  set  is  enlarged, 
numerical  convergence  of  the  results  occurs.  Furthermore,  similarities 
in  numerical  convergence  between  theoretical  approximations  establish 
criteria  for  choosing  basis  sets  in  subsequent  applications.  For  example, 
numerical  convergence  with  respect  to  basis  sets  in  the  ASCF  approxima- 
tion might  be  indicative  of  convergence  in  the  correlation  self-energy. 
For  propagator  calculations,  where  one  basis  set  describes  more  than 
one  state,  the  use  of  different  basis  sets  for  each  state  in  the  ASCF 
calculations  is  undesirable.  When  convergence  of  SCF  total  energies  for 
each  state  is  attained  with  basis  set  augmentations,  the  occupied 
orbitals  and  their  energies  have  been  adequately  characterized.  Virtual 
orbitals,  however,  may  change  markedly  as  the  basis  set  is  improved. 
Diffuse  basis  functions  have  a  pronounced  effect  on  the  virtual  orbitals. 
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The  correlation  self-energy  exhibits  the  aggregate  effect  of  these 
apparent  instabilities. 

As  the  data  in  Tables  1-2  and  1-3  demonstrate,  the  augmentations  may 
affect  the  SCF  total  energy  of  one  state  more  than  the  other.  One  p  and 
two  diffuse  s  functions  are  added  to  basis  A  to  form  basis  B.  These 
additions  cause  a  much  greater  energy  lowering  for  Cl?~  than  for  Cl?. 
Further  improvements  of  this  kind  lead  to  basis  C,  but  the  ASCF  electron 
affinities  change  little.  The  addition  of  d  functions  to  the  original 
basis  (Basis  D  and  Basis  E)  has  an  important  effect  on  the  ASCF  electron 
affinity  at  the  equilibrium  internuclear  distance  of  Cl?,  but  has  a  less 
important  effect  at  the  Cl2"  equilibrium  distance.  Results  for  the 
second  order  relaxation  contributions  to  the  self-energy  are  also  listed 
in  Table  1-4.  These  terms  change  similarly  to  ASCF  with  respect  to 
basis  set  improvements.  (See  also  the  restricted  SCF  study  of  Gilbert 
and  Wahl,  1971.) 

Second  order  and  diagonal  2p-h  TDA  correlation  self-energies  are 
constructed  with  the  orbitals  from  the  Cl?  SCF  calculation.  The  improve- 
ment of  the  correlation  description  decreases  the  vertical  electron 
affinity,  but  improved  basis  sets  increase  it. 

For  CI-,  as  with  many  molecules,  experimentalists  have  measured 
only  the  adiabatic  electron  affinity  (Dispert  and  Lacmann,  1977).  Energy 
differences,  calculated  or  empirical  (Douglas  et  al.,  1963),  along  the 
CI 2  or  CI 2  potential  curves  are  necessary.  Energy  differences  in  the 
SCF  approximation  for  the  Clp  curve,  as  well  as  second  order  RSPT  calcu- 
lations for  Cl2,are  in  Table  1-5.  Adiabatic  electron  affinities  using 
various  combinations  of  potential  energy  curves  and  vertical  electron 
affinities  are  in  Tables  1-5,  1-6  and  1-7.  The  restricted  Hartree-Fock 
potential  curve  is  inadequate  because  it  separates  into  ions  instead  of 
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Tabl 

e  1-7 

Adiabat 

ic  Electron 

Affinity 

of  CI  2 

(eV) 

Ver 

tical 
(A-1) 

v(ci2") 

EA(adiab) 

Eq- 

3.08 

1.59 

Eq. 

(A-2) 

2.62 

1.13 

Eq. 

(A-1) 

3.60 

2.11 

Eq. 

(A-2) 

3.19 

1.70 

Eq. 

(A-1) 

4.00 

2.51 

Eq- 

(A-2) 

3.74 

2.25b 

D(C1~)  =  1.49  eV  (see  Douglas  et  al . ,  1963) 


a  Eq.  (A-1)  means  that  V(C12")  is  calculated  as  -I  +  LC'2'.  Eq  (A-2) 
means  that  V(C1~~)  is  calculated  as  -I  +  £  p" 

b  EXPERIMENT  EA(adiab)  =  2.4  ±  0.2  eV  (Dispert  and  Lacmann,  1977). 
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into  neutral  chlorine  atoms.  The  second  order  RSPT  results  contain  some 
correlation.  The  experimental  energy  difference  or  the  second  order 
RSPT  energy  difference  along  the  CI-  curve,  combined  with  the  vertical 
electron  affinity  at  5.0  bohr,  give  an  adiabatic  electron  affinity  that 
is  within  experimental  error.  A  smaller  basis  set  combined  with  a 
simpler  treatment  of  correlation  could  have  produced  a  prediction  that 
was  also  within  experimental  tolerance  through  a  fortuitous  cancellation 
of  errors. 

1.7  The  electron  affinity  of  CN 

The  vertical  electron  affinity  of  CN  is  negligibly  different  from 
the  adiabatic  electron  affinity  since  there  is  little  change  in  the  bond 
length  upon  addition  of  an  electron.  (See  Table  1-11.)  The  total 
energies  listed  in  Tables  1-9  and  1-10  show  the  small  change  in  the  SCF 
total  energy  from  one  equilibrium  distance  to  the  other.  Therefore, 
vertical  electron  affinities  are  computed  at  only  one  point  per  basis 
set.  (See  Table  1-8.)  As  the  total  energies  indicate,  the  d  functions 
produce  the  greatest  improvements.  The  addition  of  diffuse  functions 
is  not  as  important  since  basis  A  has  already  been  optimized  for  anionic 
molecules  (Pacansky  and  Liu,  1977).  Griffing  and  Simons  also  studied 
the  electron  affinity  of  CN  using  the  equation  of  motion  technique 
(Griffing  and  Simons,  1976).  They  chose  Hartree-Fock  plus  first  order 
RSPT  corrections  for  the  ground  state  and  the  f-  operator  product  mani- 
fold. They  also  varied  the  exponents  of  the  s  and  p  Gaussian  basis 
functions  to  minimize  the  5a  orbital  energy  (see  the  critique  of  Liu, 
1977)  and  included  no  d  functions  whatsoever. 

The  two  smallest  basis  sets  permit  third  order  correlation  self- 
energy  calculations  (see  Table  1-12).  For  the  smallest  basis  set,  the 
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third  order  result  is  in  closest  agreement  with  experiment,  but  the 
addition  of  the  d  functions  casts  doubt  on  the  balance  of  this  approxima- 
tion. If  the  trends  of  the  basis  set  improvement  continue,  the  third 
order  approximation  will  be  inferior  to  the  diagonal  2p-h  TDA  result  as 
the  basis  set  is  improved. 

1.8  Ionization  energies  of  NH2~ 


The  lowest  ionization  energy  of  NHL"  is  the  electron  affinity  of  NH«. 
The  second  ionization  energy  of  NFL"  yields  an  excited  final  state  of  NFL 
and,  in  combination  with  the  lowest  ionization  energy,  an  excitation 
energy.  For  the  first  ionization,  there  is  little  nuclear  rearrangement 
and,  as  with  CN,  the  vertical  electron  binding  energy  at  one  nuclear 
geometry  is  sufficient.  For  the  second  ionization,  however,  there  is  a 
large  change  in  the  bond  angle  and  a  small  change  in  the  bond  distances. 
Once  again,  energy  differences  along  the  N  electron  or  N±l  electron 
surface  are  necessary.  Other  investigators  (Bender  and  Schaefer,  1971; 
Heaton  and  Cowdry,  1975)  have  calculated  energy  differences  using  SCF 
and  CI  wavefunctions.  These  calculations  indicate  that  SCF  calculations 
are  reliable  for  the  surfaces  of  interest.  (Heaton  and  Cowdry  calculated 
ASCF  electron  binding  energies  of  NFL  with  a  basis  set  optimized  only  for 
uncharged  molecules.)  The  first  basis  consists  of  the  nitrogen  5s,  4p, 
Id  basis  from  the  CN"  study  plus  the  3s,  lp  hydrogen  basis  of  Table  1-13. 
The  second  basis  consists  of  an  8s,  7p,  2d  nitrogen  basis  (Dunning,  1970) 
and  a  7s,  lp  hydrogen  basis  (Huzinaga,  1965).  (See  Table  1-14.)  The 
clear  superiority  of  the  second  basis  over  the  first  does  not  manifest 
itself  by  large  changes  in  the  electron  binding  energies.  The  differences 
for  all  theoretical  approximations  between  the  basis  sets  are  small  and 
indicate  that  the  first  basis  is  adequate.  [Anderson  and  Simons  (1976) 
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Table  1-14 
CGTO  Basis  for  NH, 

NITROGEN3 


Type 

Exponent 

Coefficient 

Type 

Exponent 

Coefficient 

s 
s 

13520. 
1999. 
440. 
120.9 
38.47 
13.46 

13.46 
4.993 

0.000760 
0.006276 
0.032847 
0.132396 
0.393261 
0.546339 

0.252036 
0.779385 

P 

P 
P 
P 

35.91 
8.48 
2.706 

0.9921 

0.3727 

0.1346 

0.040319 
0.243602 
0.805968 

1.0 

1.0 

1.0 

s 

1.569 

1.0 

P 

0.059 

1.0 

s 

0.58 

1.0 

P 

0.033 

1.0 

s 

0.1923 

1.0 

P 

0.21 

1.0 

s 

0.088 

1.0 

s 

0.048 

1.0 

d 

1.5 

1.0 

s 

0.30 

1.0 

HYDROGEN 

d 
b 

0.95 

1.0 

s 

68.16 

1.0 

s 

0.082217 

1.0 

s 

10.2465 

1.0 

s 

0.03580 

1.0 

s 

2.34648 

1.0 

s 

0.673320 

1.0 

P 

1.0 

1.0 

s 

0.22460 

1.0 

a  Dunning,  1970. 
b  Huzinaga,  1965. 
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found  that  d  functions  are  important  in  describing  the  ground  state 
potential  energy  surface,  but  their  effect  on  equation  of  motion  calcula- 
tions for  the  electron  affinity  of  NH2  is  small.]  The  first  basis  is 
small  enough  to  permit  a  third  order  correlation  self-energy  calculation. 
(See  Table  1-15.) 

The  results  of  the  calculations  indicate  that  the  diagonal  2p-h  TDA 
is  the  most  reliable  predictor.  Although  the  third  order  results  are  not 
markedly  inferior  for  the  vertical  ionization  energy  of  NH2  ,  they  give 
less  accurate  values  for  the  adiabatic  splitting  of  the  NH2  states. 
(See  Table  1-16.)  The  adiabatic  electron  affinity  obtained  with  the 
anionic  surface  is  slightly  better  than  that  estimated  with  the  neutral 
surface.  The  satisfactory  approximation  of  the  electron  affinity  of  NH2 
and  the  energy  splitting  of  the  NHL  states  indicate  that  the  ASCF  plus 
d- agonal  2p-h  TDA  correlation  self-energy  is  a  balanced  treatment  of 
electron  affinities. 

1.9  Ionization  energies  of  PH2~ 

The  states  of  PH2"  and  PH2  are  similar  to  those  of  NH2"  and  NH2< 
These  calculations  show  if  the  approach  taken  for  NHL"  with  regard  to 
the  basis  set  choice  and  to  the  treatment  of  correlation  is  valid  for  a 
molecule  with  a  third  period  atom.  The  phosphorus  basis  set  (McLean  and 
Chandler,  1980)  is  approximately  as  good  as  the  first  basis  set  for  NH2". 
(See  Table  1-17.)  The  phosphorus  basis  set  has  been  optimized  for  the 
anionic  atom  as  was  the  first  nitrogen  basis  of  the  NhL  study.  The 
hydrogen  basis  is  in  Table  1-13.  The  nuclear  geometries  of  the  neutral 
states  and  of  the  anionic  state  are  optimized  with  this  basis.  The  third 
order  self-energy  is  certainly  no  better  than  the  diagonal  2p-h  TDA 
self-energy.  The  agreement  with  experiment  for  the  vertical  electron 
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Table  1-15 


Vertical  Ionization  Energies  for  NH  ~  {e\l] 


Basis 

Bond  Distance 
R  (bohr) 

Bond  Angle 
9  (deg.) 

Symmetry 
i 

"£i 

N:  5s,4p,ld 
H:  3s, lp 

1.92 

103.9 

ibj 

3a, 

1.28 
3.63 

N:  8s,7p,2d 

1.92C 

102. 6C 

ib, 

1.32 

H:  7s, lp 

3a, 

3.87 

a  Compare  difference  with  the  observed  splitting  of  2A,  and  2B,  (1.98  to 
b  Compare  with  experiment  (0.779  ±  0.037  eV);  Celotta  et  al .  (1974). 
c  Heaton  and  Cowdry  (1975). 
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Table  1-15  -  Extended 


-(-I.j(ASCF) 

I.(ASCF) 

+^2>ii(E)) 

-1.10 

0.21 

1.03 

2.22 

-1.09 

0.31 

1.11 

2.36 

-(-^•(ASCF)  -(-I.(ASCF: 

+zC(2p-hTDA).,(E))     +zC(3),n.(E); 


2.4  eV);  Dressier  and  Ramsay  (1959). 
Anderson  and  Simons  (1976)  obtain  0.42  eV. 


n 


0.61b  0.52 


2.50  3.56 

0.67a'b 
2.63a 


1—1         I/) 


S_  i—i 
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Table 

1-17 

CGTO  Basis 

for 

PH2" 

PHOSPHORUS9 

Type 

Exponent 

Coefficient 

Type 

Exponent 

Coefficient 

5 

77492.4 

0.00781 

P 

384.84 

0.009206 

11605.8 

0.006068 

90.552 

0.069874 

2645.96 

0.031160 

29.134 

0.292470 

754.98 

0.123431 

10.886 

0.728103 

248.76 

0.378209 

91.157 

0.563262 

P 

4.3526 
1.7771 

0.628349 
0.428044 

S 

91.157 

0.160255 

36.226 

0.627647 

P 

0.6970 

1.0 

15.211 

0.263849 

P 

0.2535 

1.0 

s 

4.7942 

1.0 

P 

0.0685 

1.0 

s 

1.8079 

1.0 

s 

0.3568 

1.0 

d 

0.43C 

1.0 

s 

0.1148 

1.0 

HYDROGEN 
see  Table  1-131 


a  McLean  and  Chandler,  1980 

b  Huzinaga,  1965. 

c  Dunning  and  Hay,  1977. 
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affinity  of  PH2  is  not  as  good  as  with  NH2  (see  Table  1-18).  The 
adiabatic  splitting  of  the  PH2  states  is  satisfactory  in  the  diagonal 
2p-h  TDA,  although  the  third  order  approximation  is  not  so  far  off  this 
time  as  it  was  with  NFL.  (See  Table  1-16.) 

The  computational  experience  gained  in  these  studies  indicates  that 
ASCF  +  diagonal  2p-h  TDA  self-energy  is  a  useful  tool  for  studying 
electron  affinities.  It  is  sufficiently  simple  to  allow  application  to 
larger  molecules. 

1.10  The  transition  operator  method 

For  considering  the  electron  binding  energies  corresponding  to  the 
attachment  or  removal  of  an  electron  in  spin-orbital  i,  a  choice  of  spin- 
orbital  occupation  numbers  (Purvis  and  Ohrn,  1976), 

(  0  ,  j  unoccupied 

nj  =  1*5  .  J  =  i  (1-44) 

[  1  ,  j  occupied 

defines  a  grand-canonical  one-electron  reduced  density  matrix  (Abdulnur 
et  al . ,  1972).  This  density  matrix  gives  the  form  of  a  transition 
operator  matrix  with  the  elements 

Tij  =  hij  +  I   <ik||jk><nk>  (1-45) 

where  the  subscripts  are  spin-orbital  indices.  This  matrix  is  set  up  and 
diagonal ized  iteratively  in  much  the  same  manner  as  the  Fock  matrix  in 
UHF  spin-orbital  optimizations.  At  convergence,  the  spin-orbitals 
satisfy 
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The  i-th  spin-orbital   energy  agrees  with  ASCF  through  third  order  except 
for  some  correction  terms   (Goscinski   et  al.,   1973;   Born  et  a!.,   1978). 
These  terms  have  been  evaluated  in  a  number  of  systems  and  are  small. 
(They  are  less  than  0.01  electron  volts  for  all   the  basis  sets  of  CN~ 
and  Cl?.)     This  method  incorporates  many  relaxation  effects  into  a  single 
spin-orbital  and  suggests  a  reference  state  for  a  propagator  calculation. 
Specifically,  the  transition  operator  method  (TOM)  spin-orbitals  and 
energies  are  the  basis  of  the  self-energy.     The  ground  state  density  is 
the  grand-canonical   density  formed  from  the  occupation  numbers  of 
equation  1-44.     It  is  possible,  then,  to  evaluate  second  order  and 
diagonal   2p-h  TDA  corrections  to  the  TOM  spin-orbital   energy  in  the 
quasi-particle  approximation,  but  without  excluding  terms  in  the  self- 
energy  that  correspond  to  relaxation  when  Hartree-Fock  orbitals  are 
used.     Since  the  relaxation  is  itself  renormalized  into  the  orbital, 
there  is  no  need  to  restrict  the  summations  in  the  self-energy. 

1.11     Applications  of  the  transition  operator  reference 

TOM  spin-orbital   energies  and  spin-orbital   energies  plus  self-energy 
corrections  are  in  Table  1-20,  along  with  results  in  the  same  basis  for 
H20  (Dunning,   1970;  Huzinaga,   1965).     (See  Table  1-19.)     Among  the  other 
approximations  are  the  ASCF  plus  correlation  self-energy  results.     The 
TOM  plus  diagonal   2p-h  TDA  self-energy  results  are  superior  to  the  other 
approximations. 

For  Ne,  the  TOM  plus  self-energy  calculations  employ  a  large  Gaussian 
basis  set  (Table  1-21).     While  the  2s  and  2p  ionizations  are  satisfacto- 
rily predicted,  the  Is  prediction  is  much  poorer.     The  diagonal   2p-h  TDA 
result  is  much  worse  than  the  second  order  result  (see  Table  1-22);   this 
effect  could  be  due  to  the  necessity  for  precisely  describing  certain 
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Table  1-19 
26  CGTO  Basis  for  H20 

Type  Exponent  Coefficient 

0.032828 

0.23121 

0.81724 

1.0 

1.0 

0.002031 

0.015436 

0.073771 

0.2476 

0.61183 

0.2412 

1.0 

1.0 

1.0 

0.01958 

0.12419 
0.39473 
0.6273 

1.0 

1.0 


H,s 

17.37 

2.6273 

0.58994 

H,s 

0.16029 

H,p 

1.0 

0,s 

7816.5 

1175.8 

273.19 

81.17 
27.184 
3.4136 

0,s 

9.5322 

0,s 

0.9398 

0,s 

0.2846 

0,p 

35.183 

7.904 

2.3051 

0.7171 

0,p 

0.2137 

0,d 

1.0 

Oxygen  s  and  p  basis  from  Dunning,  1970. 
Hydrogen  s  basis  from  Huzinaga  (1965). 
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Table  1-20 

Vertical  Ionization  Energies  of  FLO  (eV] 
26  CGTO  Basis  Set 


lal       2al      lb2      3al      lbl 


TOM 

540.35 

33.97 

17.41 

13.03 

10.75 

TOM  +  2nd  order 

540.51 

33.36 

18.56 

14.13 

11.84 

TOM  +  diag.  2p-h 

TDA 

33.56 

18.94 

14.61 

12.38 

Orbital  energy 

559.39 

36.62 

19.35 

15.66 

13.67 

ASCF 

540.49 

33.86 

17.93 

13.15 

10.89 

ASCF  +  2nd  order 

540.70- 

32.78 

18.99 

14.20 

11.93 

ASCF  +  diag.  2p-r 

TDA 

540.62 

31.28 

19.20 

14.44 

12.19 

Second  order 

33.4 

18.0 

13.4 

11.1 

Diag.  2p-h  TDAa 

33.7 

18.5 

14.1 

11.8 

Second  order 

32.93 

17.70 

13.18 

10.92 

Third  order 

35.10 

19.22 

15.18 

13.03 

Experiment 

540.2 

32.2 

18.6 

14.7 

12.6 

a  Diagonal ization  of  full  propagator  matrix  (Born  and  Ohrn,  1979). 
b  Diagonal ization  of  full  propagator  matrix  (Cederbaum  et  al.,  1973) 
c  Photoelectron  spectroscopy  (Siegbahn  et  al. ,  1969) . 
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Table  1-21 
Neon  Gaussian  Basis  Set3 

Type  Exponent 

s  279341.011322 

43638.324426 
9155.291348 
2399.53 
744.787 
260.221 
100.115 
41.7663 
18.6344 
8.76275 
4.21282 
1.7713 
0.723060 
0.285620 

p  354.616931 

82.0552 
25.6929 
9.44944 
3.82817 
0.70396 
0.30883 
0.13229 

d  2.5 

2.0 
0.8 


a  The  basis  employs  no  contractions. 

See  van  Duijneveldt  (1971)  for  details 
of  basis  set  optimization. 
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Table  1-22 
Ionization  Energies  for  Ne  (eV) 


Is.  2s_  2£ 


TOM 

868.14 

49.09 

19.54 

TOM  +  2nd  order 

868.98 

48.37 

20.91 

TOM  +  diag.  2p-h  TDA 

873.1 

48.57 

21.58 

Second  order 

19.86 

Third  order 

22.40 

Valence  approximation 

21.56 

Experiment 

870.2 

48.26 

21.56 

a  Photoelectron  spectroscopy  (Carlson  et  al . ,  1971) 
b  Cederbaum  and  Domcke  (1977). 


45 


Table  1-23 


Vertical  Ionization  Energies  of  NH?~  and  CN~ 

Calculated  with  the  Transition  Operator  Method 

and  Self-Energy  Corrections  (eV) 


System    Orbital     TOM     TOM  +  2nd  order  TOM  +  diaq  2p-h  TDA 

2.80  3.57              3.86 

0.78  1.93              2.51 

-1.37  -0.24              0.30 

CN":  Basis  D  from  Table  1-8. 
NH2":  First  Basis  from  Table  1-15. 


CN" 

5a 

NH2" 

3al 

NH2" 

lbl 

46 
integrals  involving  core  orbitals  which  arise  in  the  denominator  shift. 
Such  precision  may  be  lacking  when  a  Gaussian  basis  set  is  employed  (see 
Cederbaum  and  Domcke,  1977). 

In  the  case  of  the  electron  affinity  of  CN,  the  results  using  the 
largest  basis  set  from  the  previous  group  of  calculations  are  superior 
to  all  other  results.  (See  Table  1-23.) 

Finally,  less  satisfactory  results  are  obtained  for  the  NH2"  ioniza- 
tions. Since  ASCF  calculations  are  worse  than  Koopmans's  theorem  predic- 
tions in  this  case,  the  relaxation  renormalization  may  hinder  the 
convergence  of  self -energy  improvements. 

1.12  Interpretation 

Quasi-particle  approximations  are  conservative  methods  for  extending 
the  orbital  energy  concept  to  explain  electron  binding  energies  without 
invoking  the  full  apparatus  of  electron  propagator  theory.  The  overlap 
amplitudes  remain  proportional  to  simple  Hartree-Fock  canonical  orbitals. 
There  are  many  poles  that  correspond  to  each  orbital,  but  the  poles  have 
different  weights.  Thus,  certain  so-called  shakeup  peaks  in  a  photoelec- 
tron  spectrum  may  be  assigned  to  a  simple  orbital.  The  success  of  these 
methods  implies  that  it  may  not  be  necessary  to  completely  abandon  the 
simple  orbital  picture  to  explain  binding  energies,  particularly  in  the 
core  and  outer  valence  regions. 


CHAPTER  TWO 
EXCITATION  ENERGIES 

2.1  Consistent  ground  states  and  the  random  phase  approximation 

The  previous  chapter's  study  of  the  uncorrected  electron  propagator 
and  its  extension  to  correlated  electron  dynamics  is  an  example  of  how 
approximate  propagators  are  constructed  to  give  improved  energy  differ- 
ences and  transition  probabilities.  Avoiding  the  construction  of  wave- 
functions  for  states  other  than  the  ground  state  is  the  salient  differ- 
ence between  the  propagator  approach  to  correlation  and  variational 
theory.  This  characteristic  is  at  once  the  strength  and  the  enigma  of 
propagator  theory  (Linderberg  and  Ohrn,  1977).  Since  improved  ground 
state  density  matrices  from  the  propagator  are  useful  for  the  calculation 
of  properties  (Jorgensen,  1975;  Purvis  and  Ohrn,  1975),  it  is  important 
to  know  if  there  exists  a  corresponding  wavefunction.  This  is  the  so- 
called  N-representability  problem  (Davidson,  1976  and  references  therein) 
which  also  occurs  in  other  theories  that  seek  to  circumvent  the  explicit 
use  of  wavefunctions,  such  as  density  matrix  and  density  functional 
theories  (Ohrn,  1976).  The  question  of  N-representability  in  propagator 
theory  can  be  formulated  in  terms  of  so-called  consistency  requirements 
or  in  terms  of  sum  rules,  which,  when  satisfied,  imply  the  existence  of 
wavefunctions  for  the  ground  and  excited  states  generated  by  the  propaga- 
tor (Linderberg  and  Ohrn,  1968). 

To  elucidate  the  points  made  above,  the  model  electron  propagator 
of  the  previous  chapter  will  serve  as  an  example  of  a  consistent 
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propagator.  The  model  electron  propagator  with  the  Hartree-Fock  ground 
state  and  the  operator  manifold  consisting  of  simple  field  operators  has 
electron  addition  operators 

|N+1>  =  a  +|N>  (2-1) 

which  generate  N+l  electron  states  and  electron  annihilation  operators 

|N-1>  =  ah|N>  (2-2) 

which  generate  N-l  electron  states.  The  adjoints  of  these  operators 
annihilate  the  ground  state,  i.e., 


ap|N>  =  ah+|N>  =  0  .  (2-3) 


Thus,  the  model  electron  propagator  provides  an  explicit  wavef unction  for 
the  ground  state  and  the  N+l  electron  states.  If  the  model  electron 
propagator  is  constructed  in  a  basis  so  that  its  inverse  matrix  elements 
are 


VE)  =  E5ij  "  (hij  +  j[I<iM|Jl><akV  -  (El-F)^  ,   (2- 


4) 


(see  equations  1-5  for  an  explanation  of  the  symbols),  then  this  matrix 
is  easily  related  to  the  Fock  matrix  which  is  constructed  in  a  Hartree- 
Fock  self-consistent  field  calculation  (see  McWeeny  and  Sutcliffe,  19691 
When  the  eigenvectors  of  this  propagator  are  reinserted  into  the 
propagator,  an  identical  self-consistent  calculation  ensues.  The 
Hartree-Fock  SCF  procedure  is  a  self-consistent  electron  propagator 
calculation.  This  self-consistent  procedure  is  justified  because  it  is 
known  that  the  consistent  ground  state  for  the  simple  field  operator 
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manifold  is  a  determinantal  wavefunction.  Consistent  ground  states  are 
states  where  the  adjoints  of  the  excitation  operators  (or  ionization 
operators  in  the  case  of  the  electron  propagator)  destroy  the  ground 
state.  If  the  operator  manifold  is  extended  beyond  simple  field  opera- 
tors, then  the  consistency  requirements  are  no  longer  fulfilled.  The 
Dyson  equation  will  yield  correlated  binding  energies  and  overlap 
amplitudes,  but  there  will  be  no  explicit  wavef unctions  for  the  N±l 
electron  states  and  no  consistent  ground  state.  The  propagator  can  yield 
an  improved  ground  state  one-electron  density  matrix,  but  it  is  not 
necessarily  N-representable. 

A  correlated  consistent  ground  state  is  associated  with  the  particle- 
hole  or  polarization  propagator.  This  propagator  produces  excitation 
energies  and  overlap  amplitudes  between  N  electron  states  as  the  electron 
propagator  produces  ionization  energies,  electron  affinities  and  overlap 
amplitudes  between  the  ground  state  and  N±l  electron  states  (Linderberg 
and  Ohrn,  1973).  The  polarization  propagator,  expressed  as  a  matrix  in 
a  basis  of  spin-orbitals,  is  given  by  the  expression 


pijki<E= 


i  t  i    1  +  1 
<0 1  a .  a, |n><n|a.  a, |0> 


E-En+E0+ic 


4.  4. 

<0 1  a.  a, |n><n|a.  a . |0> 


E+En"E0- 


la 


!2-5) 


Note  that  the  poles 


E  =  *<En-E0> 


:2-6! 


and  the  residues 


4-  4. 

<0|a.  a  .|nxn|a.  a-,  |0> 


2-l\ 


have  an  interpretation  for  excitation  processes  that  is  completely 
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analogous  to  the  interpretation  of  the  poles  and  residues  of  the  electron 
propagator  for  ionization  processes.  As  the  electron  propagator  was 
developed  to  give  the  Dyson  equation  from  the  inner  projection  technique, 
the  polarization  propagator  can  yield  the  Bethe-Salpeter  equation  with  the 
aid  of  inner  projection  theory.  The  simplest  operator  manifold  that  is 
inspired  by  the  inner  projection  technique  uses  the  f~  operator  product 
manifold: 


P(E)  =  (a+a|h)(h|(EI-H)h)"1(h|a+a)  ,  (2-8) 


where 


h  =  f2  =  {ak+ark/l}  .  (2-9) 

Improvements  can  be  realized  by  including  higher  operator  products,  such 
as  f^  or  fg.  Now  the  question  of  the  ground  state  choice  arises.  The 
choice  of  the  Hartree-Fock  ground  state  produces  the  time-dependent 
Hartree-Fock  approximation  (Jorgensen  and  Linderberg,  1970).  This  ground 
state  is  not  consistent  and  leads  to  inquiries  concerning  the  form  of  the 
ground  state  that  is  consistent  with  f2  (Linderberg  and  Ohrn,  1977;  Ohrn 
and  Linderberg,  1979;  Goscinski  and  Weiner,  1980;  Weiner  and  Goscinski, 
1980a  and  1980b).  The  form  of  the  consistent  ground  state  is  the  anti- 
symmetrized  geminal  power  (AGP), 

N/2  , 
|ground  state>  =  0»<-  h  g(2i-l,2i)  ,  (2-10) 

rti  i=l 

where  N  is  the  number  of  electrons  and  0AS  =  jjjy  I   (-1)PP,  which  suggests 
that  a  variationally  optimized  AGP  might  be  an  interesting  ground  state 
for  a  polarization  propagator  calculation  with  the  f~  operator  product 
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manifold.  All  polarization  propagator  calculations  that  employ  the  f? 
operator  manifold  are  called  random  phase  approximation  calculations. 
This  name  is  also  applied  when  a  Hartree-Fock  reference  state  is  used. 
The  name  extended  RPA  is  applied  to  calculations  where  there  are  improve- 
ments in  the  ground  state  or  in  the  operator  manifold  (Shibuya  and  McKoy, 
1970).  Finally,  the  term  self-consistent  RPA  has  been  applied  to  calcu- 
lations where  improved  ground  state  densities  from  the  polarization 
propagator  have  been  calculated  iteratively  (Jorgensen  and  Oddershedde, 
1972).  In  the  present  calculations  there  is  no  assurance  that  the 
variational^  determined  AGP  will  be  the  ground  state  that  fulfills  all 
the  consistency  requirements.  Therefore,  the  polarization  propagator 
calculations  that  follow  are  called  AGP-RPA. 

The  AGP  has  previously  been  studied  for  reasons  that  are  unrelated 
to  propagator  theory  (Kutzelnigg,  1964  and  1965;  Bratoz  and  Durand,  1965; 
Bessis  et  al . ,  1967  and  1969;  Coleman,  1965;  Fukutome,  1977;  Fukutome  et 
al . ,  1977).  The  sparseness  of  its  two-electron  density  matrix  and  the 
incorporation  of  correlation  into  a  single  geminal  have  made  AGP  attrac- 
tive to  many  investigators.  The  following  discussion  treats  the  mathe- 
matical properties  of  the  AGP  and  the  interpretation  of  electronic 
structure  that  it  provides.  In  addition,  the  AGP  is  used  as  a  reference 
state  in  a  polarization  propagator  calculation. 

2.2  The  anti symmetrized  geminal  power  wavefunction 
The  AGP 

N/2 
Cv  n  g(2i-l,2i)  (2-11) 

rt:>  i=l 

is  formed  from  the  geminal,  g(l,2),  which  is  expressed  in  terms  of  a 
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spin-orbital  basis  as  follows, 

M   M 


9(1-2)  =  I       I     C..[$.(l)0'.(2)  -  $',(1)4.(2)]  (2-12; 

1=1  j=l  1J   ]     J       J    i 

where  M  is  the  dimension  of  the  orbital  basis.  This  geminal  can  also  be 
written  in  natural  form 


M 
9(1.2)  =   J  g,[«,(D*'i(2)  -  *',(1)$,(2)]  (2-13; 

1=1         1  1    ' 


where  unprimed  $  stands  for  a  spin-orbital  with  alpha  spin  and  primed  $ 
stands  for  a  spin-orbital  with  beta  spin.  The  Hartree-Fock  wavefunction 
is  included  as  a  special  case  of  the  AGP  by  restricting  M  to  be  the 
number  of  electrons  divided  by  two.  For  the  case  of  odd-electron 
systems,  the  generalized  AGP 

N+K       N/2 
0Ac   n   Xi(i)  n  g(2i-l,2i)  .  (2-14) 

Ab  i=N+l  n    i=l 

where  the  system  has  N  +  K  electrons,  is  the  wavefunction  that  fulfills 
the  consistency  requirements  of  the  RPA.  Returning  to  the  simple  AGP, 
a  superposition  of  configurations  is  generated  when  M  is  greater  than 
half  the  number  of  electrons.  This  expansion  of  configurations  is 


I    I     ...  I 

W       >iN/29il9i2'"9iN/2     X 

det{$.    (1)$'      (2)*.    (3)*'.    (4)...$,        (N-l)*'  (N)}  (2-15' 

1  1         n2  32  \/2  nN/2 

and  includes  certain  double,  quadruple,  hextuple  and  higher  even-fold 
excitations  from  any  determinant  which  is  chosen  as  the  reference 
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configuration.  In  this  study,  the  geminal  will  be  optimized  in  its 
natural  form  and  the  further  restrictions  to  the  singlet,  restricted  case 


M  . 

9(1,2)  =  I     9,^(1)^(2)  1(a(D3(2)-6(l)a(2))  (2-16; 

i=l  /2 


will  be  made,  where  {*•}  stand  for  spatial  orbitals. 


2.3  Density  matrices  and  the  total  energy  expression 
The  singlet,  restricted  geminal, 


M   M  . 

9(1.2)  =  I       I     Ci.4.(l)*-(2)-^(a(l)6(2)-B(l)a(2))  ,      (2-17; 
i=l  j=l  1J  1    J   /2~ 


may  be  transformed  to  natural  form, 


M  . 

9(1,2)  =  I     g.<b.(lU.(2)-±(cx(l)e(2)-3(l)a(2))  ,  (2-18] 

i=l         JZ 


where  the  natural  orbitals  of  the  geminal  are  {<)>.}  and  the  natural 

2 

orbital  occupation  numbers  of  the  geminal  are  |g.|  .  The  superposition 

of  configurations  derived  from  the  AGP  contains  no  configuration  that  is 
singly  excited  with  respect  to  any  other  determinant  when  the  geminal  is 
in  natural  form.  The  elements  of  the  one-electron  density  matrix  of  the 
AGP 


N..  -  <aitaj>  -  Nii6ij  (2-19) 


are  zero  unless  the  indices  are  equal.  Therefore,  the  natural  orbitals 
of  the  AGP  are  the  natural  orbitals  of  the  geminal.  The  diagonal  element 
N . .  is  the  sum  of  the  squared  absolute  values  of  the  coefficients  of  the 
configurations  in  which  the  orbital  i  appears: 
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"ii  =  l9i|2  11      ...  I     l9i  [  2i  3-i  |2...|9i    I2  •      (2-20) 
VV   >1N  .    !    2      >-l 


By  introducing  the  symmetric  function  of  order  N/2  of  the  squares  of  the 
coefficients,  g.,  the  occupation  numbers  of  the  one-electron  density 
matrix  of  the  geminal, 

C[f)  -     I     I      •••!   191  l^i  |2---l9i   I'  •  (2-21) 

U     11>12>  \/2       l  2       N/2 

the  elements  of  the  one-electron  density  matrix  of  the  AGP  can  be 
expressed  as 


1     9C \?\ 
I..   =  — —  q.  — LfJ-  .  (?-??) 


The  elements  of  the  two-electron  density  matrix  that  are  not  zero  follow 
in  a  similar  way  according  to  the  formulas 


t     t                         1                9  C  2 
<a.   a.   a-a.>     =    — —  q.q.  L£-L  f?-?V 


C!2 


where  i  and  j  are  the  spin-orbital s,  and 


Ml+i 


t  t         1   °  °  7 

<ai  a!  VJ*  =  lo  ~fe"  '  (2-24: 


where  i  =  i  with  opposite  spin.  The  total  energy  is 
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E  =  I   (2h,.+<ii|ii>)N..  +  J  I   (4<i j | ij>-2<i j | ji>)<a,Ta,Ta  .a.>    (2-25; 
i  i  >  j  |  j  j 

+  2Re{J  I   <ii|jj><a,Ta,Ta,a,>} 
i>j  -  *  J 

(Linderberg,  1980),  with  i  and  j  now  spatial  orbital  indices. 

2.4  Interpretations  of  the  total  energy  expression 

The  total  energy  expression  can  be  interpreted  with  the  terminology 
of  configuration  interaction.  The  first  two  summations  in  equation  2-25 
arise  from  Hamiltonian  matrix  elements  between  identical  configurations. 
Each  configuration's  total  energy  is  weighted  by  the  absolute  value 
squared  of  its  coefficient  in  the  CI  expansion  generated  by  the  AGP. 
Thus,  the  integrals  that  occur  in  the  first  two  terms  are  one-electron, 
Coulomb  and  exchange  integrals.  In  the  third  summation,  the  integrals 
cannot  be  classified  as  Coulomb  or  exchange  integrals;  they  arise  from 
matrix  elements  in  the  CI  expansion  between  configurations  that  differ 
by  two  spin  orbital s. 

Another  framework  in  which  the  total  energy  may  be  analyzed  is  in 
terms  of  pair-correlation  theories  (Meyer,  1977  and  references  therein) 
and  many-body  perturbation  theory  (MBPT),  in  which  pair  correlation 
energies  provide  corrections  to  the  Hartree-Fock  total  energy  (Bartlett 
and  Shavitt,  1977;  Bartlett  and  Purvis,  1978).  These  pair  correlation 
terms  are  produced  by  matrix  elements  of  double  excitations  with  the 
closed-shell  Hartree-Fock  configuration.  In  MBPT,  the  contributions 
from  the  first  and  second  summations  in  equation  2-25  that  do  not  come 
from  the  reference  configuration  vanish  according  to  the  linked-cluster 
theorem  (Brueckner,  1955a  and  1955b;  Goldstone,  1957).  The  contributions 
from  the  third  term  derive  from  the  same  type  of  CI  matrix  elements  as 
the  pair  energies  of  the  various  many-body  approaches. 
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2.5  Coefficient  optimization 

Because  of  the  simplifications  in  the  two-electron  density  matrix 
that  obtain  by  considering  the  natural  form  of  the  geminal,  both  the 
coefficients,  g.,  of  the  geminal  and  the  natural  orbital s  of  the  geminal 
must  be  optimized.  If  the  general  form  of  the  geminal  had  been  employed 
(see  equation  2-17),  it  would  have  been  sufficient  to  only  vary  the 
coefficients,  c.  .. 

The  coefficient  variation  procedure  is  obtained  by  imposing  a 
necessary  condition  for  the  stationarity  of  the  total  energy.  Let  the 
Rayleigh-Ritz  variational  quotient, 


F  =  <AGP[  H[  AGP>  (9  ?^ 

L     <AGP|AGP>  '  {        °' 


be  written  as 


E  =  -£-   ,  (2-27; 

4 


where  E'  is  the  unnormalized  AGP  total  energy,  and  C(N/2)  is  the  normal- 
ization, or  the  symmetric  function  of  order  N/2  of  the  squares  of  the 
absolute  values  of  the  coefficients  of  the  geminal.  If  the  first  varia- 
tion of  the  energy  with  respect  to  the  coefficients  is  zero, 

E  =  0  ,  (2-28; 


3g*. 


subject  to  the  constraint  that  the  geminal  be  normalized, 


I   g*j9j  =  1  .  (2-29) 


then  the  necessary  condition 
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-  [E  "  MI  9*^-1)]  =  0 
i      j   J  J 


(2-30; 


is  obtained,  where  A  is  a  Lagrange  multiplier.  By  noticing  that  both  the 
normalization  of  equation  2-27  and  the  unnormalized  total  energy  satisfy 
the  Euler  formula  for  symmetric  functions, 


1   r  .   3E' 

n/2  4  9i  agi 


eft 


NT?  I   9- 
J 


acf 


(2-31) 


>32) 


the  necessary  condition  for  stationarity  of  the  total  energy  can  be 
rewritten 


1  3E1 


3C 


3g*-i 


oft 


3g' 


1    B# 


Nh2     y   i 


9g*,- 


(2-33) 


M?P 


2  32C&]    1 

3^E'        _   E   3  Ll2 


3g*i3gj        9g*i3gj 


b  ?  Hij9i  . 


W!"  J 


(2-34) 


(2-35) 


where  H-.  -   H*-..  Now  with  the  introduction  of  the  constraint  on  the 
normalization  of  the  wavefunction,  equation  2-30  implies  that 


f  Ml   -i    ^  J   v. 


WfN 


x9i 


!2-36; 


Equation  2-36  is  a  standard  eigenvalue  problem  that  is  iteratively 
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solved.     At  self-consistency,  the  value  of  A  is  zero  since 


"TV  2  HiigT  =  3p~  =    °  '  {2"37) 


provided  that  g.  is  itself  not  zero. 

2.6  Orbital  optimization 

A  unitary  transformation  of  the  orbital s  in  the  AGP  may  be  written  as 

ar  =  e   ar  e  (2-38) 

(Dalgaard  and  Jorgensen,   1978),  where 

A     =     i  Ipkl(ak+ai   -ai+ak)    .  (2-39) 

The  resulting  total  energy  after  such  a  transformation  on  |0>  is 

E(A)  =  <0|e"lA  HelA|0>  =  <0|H|0>  +  -i<0| [A,H] |0>         (2-40) 

+  %{<0|[A,[H,A]]|0> 
+  <0|[[A,H],A]|0>}  +  ...  . 

A  stationary  point  occurs  when  <5E(A)  =  0,  or 

6E(X)  =  -i<0|[6A,H]|0>  +  %{<0|[6A,[H,A]]|0>  +  (2-41) 

+  <0|[[6A,H],A]|0>}  +  ...  . 


Let 


59 
Wrs  =  <0|[artas  -  as+ar,H]|0>  (2-42) 


and  let 


ArSjr,s,  -  ^<0|[ar+as  -  as+ar,[H,ar,+as,  -as,+ar,]|0>    (2-43; 
+  <0|[[ar+as  -  as+ar,H],ar,+asl  -  as,+ar,]|0>  . 


With  the  neglect  of  terms  higher  than  second  order,  equation  2-41  becomes 


X     =     -A_1W  .  (2-44) 


The  solutions,  A,  are  increment  vectors  and  are  converted  to  unitary 
orbital  transformation  matrices  according  to 


X  =  e"A  =  -Ucosh(d)U+  -  AUsinh(d)d-1U+  (2-45) 


where 


A2  =  UdU+  .  (2-46) 


The  integrals  are  transformed  to  the  new  basis  before  the  next  orbital 
or  coefficient  iteration. 

2.7  Results  for  lithium  hydride 

Calculations  for  AGP  wavefunction  optimizations  are  attempted  with 
a  small  Gaussian  basis  set  (see  Table  2-1).  The  total  energies  from 
these  calculations  and  the  total  energy  differences  along  the  ground 
state  potential  energy  curve  are  compared  to  results  from  Hartree-Fock, 
second  order  MBPT,  an  MC-SCF  wavefunction  that  contains  the  same  config- 
urations as  AGP,  but  not  the  constraints  on  the  coefficients,  and  to  the 
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Table  2-1 
LiH  Basis 

13  CGTO  Basis' 


i/Ri 


Exponent 


Coefficient 


Li,s 


642.41895 

96.79849 

22.09109 

6.20107 

1.93512 

0.63674 


0.00214 
0.01621 
0.07732 
0.24579 
0.47019 
0.34547 


Li,s 


2.19146 
0.59613 


0.03509 
0.19123 


Li  ,s 
Li  ,s 


0.07455 

0.02079 
0.00676 


1.00000 

0.39951 
0.70012 


Li,s 
Li,p 


0.08948 

2.19146 
0.59613 
0.07455 


1.00000 

0.00894 
0.14101 
0.94535 


Li  ,p 


0.08948 
0.02079 
0.00676 


0.15559 
0.60768 
0.39196 


H,s 


18.73110 
2.82539 
0.64012 


0.00349 
0.83473 
0.81376 


H,s 


0.16128 


1.00000 


a  Kurtz  and  Ohrn,  1978. 
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full  CI  result.  The  percentage  of  the  correlation  energy  is  also 
monitored. 

The  usual  starting  point  for  beginning  the  optimizations  is  the 
Hartree-Fock  wavefunction  in  the  canonical  orbital  basis.  The  choice  of 
coefficients  is  zero  or  one  according  to  whether  the  orbital  in  question 
is  occupied  or  not  in  the  Hartree-Fock  wavefunction.  (The  coefficients 
are  assumed  to  be  real.)  Coefficient  optimizations  are  begun;  usually 
five  coefficient  iterations  are  sufficient  to  begin  the  orbital  optimiza- 
tion. Only  one  orbital  optimization  is  made  before  the  next  round  of 
coefficient  optimizations  is  started  (see  Table  2-2).  The  motivation 
for  favoring  the  coefficient  optimizations  stems  from  the  comparatively 
smaller  effect  of  the  orbital  optimizations  on  the  total  energy  and  from 
the  necessity  of  performing  an  integral  transformation  after  each  orbital 
optimization.  For  the  regions  of  the  potential  energy  curve  where  the 
Hartree-Fock  determinant  dominates  the  full  CI  expansion,  this  procedure 
is  satisfactory  and  converges  in  about  20  cycles,  where  each  cycle  is 
composed  of  five  coefficient  iterations  and  one  orbital  iteration.  The 

convergence  behavior  for  other  regions  of  the  curve  grows  worse  as  the 

2     2 
doubly  excited  determinant  cu  — ►  cu  becomes  more  important  in  the  CI 

expansion  of  the  wavefunction.  The  AGP  thus  begun  also  has  a  tendency 

to  converge  to  an  excited  state.  For  this  region  of  the  potential  curve 

(beyond  R  =  4.0  bohr),  the  natural  orbitals  of  the  full  CI  result  were 

used,  and  the  AGP  calculation  converged  as  quickly  as  the  previous 

results.  The  full  CI  natural  orbitals  were  also  used  to  start  AGP 

calculations  in  the  interior  of  the  potential  curve  (within  R  =  4.0  bohr) 

as  a  check  on  the  original  results.  The  AGP  iterations  converged  to  the 

original  result.  A  typical  example  of  the  convergence  behavior  is  given 

in  Table  2-2.  The  use  of  a  damping  factor  of  0.5  in  the  coefficient 
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Table  2-2 


Convergence  Behavior  for  AGP  Wavef unction 
Optimization  for  LiH  Z     at  R  =  3.015  bohr 


Cycle3  Energy  (a.u, 


0  -7.96216 

1  -7.97133 

2  -7.97416 

3  -7.97512 

4  -7.97514 

5  -7.97515 

10  -7.97515  53658 

15  -7.97515  53664 


a  One  cycle  is  five  coefficient  iterations  followed 
by  one  orbital  iteration. 
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optimizations  has  been  found  to  assure  convergence,  although  it  is  not 
necessary  to  impose  a  damping  factor  in  all  cases.  The  convergence 
behavior  listed  in  the  table  is  probably  not  the  most  rapid  possible  with 
this  method.  The  convergence  criterion  for  the  total  energy  of  10    is 
generally  the  first  criterion  to  be  satisfied.  The  sum  of  the  squares 
of  the  orbital  transformation  matrix  elements  and  the  sum  of  the  squares 
of  the  elements  of  the  orbital  gradient  matrix, 

H   |X  |2        and        I  I   |<0|[a  +as,H]|0>|2  ,    (2-47) 
r  s   rs  r  s 

are  also  employed  as  convergence  criteria.  When  these  tolerances  are 
both  set  equal  to  10"  ,  the  former  criterion  is  invariably  the  last  to 
be  satisfied,  often  requiring  five  iterations  after  the  satisfaction  of 
the  total  energy  requirement. 

Inspection  of  the  coefficients  of  the  geminal  at  each  internuclear 
distance  reveals  a  trend  (see  Table  2-3).  Near  the  equilibrium  inter- 
nuclear distance,  two  of  the  coefficients  dominate  the  geminal.  This 
dominance  is  revealed  by  the  table  in  which  natural  occupation  numbers 
for  the  AGP  are  listed.  The  other  terms  in  the  geminal  are  quite  small 
and  have  the  opposite  sign  as  the  dominant  coefficients.  As  the  bond 
length  is  increased,  a  third  coefficient  becomes  important.  This  coef- 
ficient has  a  negative  sign  and  approximates  the  lower  of  the  two  other 
important  coefficients.  The  absolute  values  of  the  second  and  third 
largest  coefficients  decrease  as  the  distance  between  the  nuclei 
increases.  The  largest  coefficient  attains  the  maximum  value  of  one. 
The  AGP  coefficients  are  forcing  the  AGP  to  approximate  a  two  configura- 
tion wavefunction  where  both  electrons  of  the  highest  occupied  sigma 
orbital  in  the  Hartree-Fock  wavefunction  are  promoted  to  the  lowest 
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unoccupied  sigma  orbital,  i.e.,  o^  a2  —  c^2^2.  This  type  of  excitation 
is  included  in  the  AGP  wavefunction.  The  next  most  important  configura- 
tion in  the  AGP  is  the  configuration  where  all  four  electrons  are 
assigned  to  the  sigma  orbitals  that  were  just  mentioned:  o^o-} .     The 
core  orbitals  on  lithium  are  not  occupied  in  this  configuration.  In  the 
limit  of  infinite  separation  of  nuclei,  such  a  configuration  should  have 
a  coefficient  of  zero  in  the  CI  wavefunction  since  it  corresponds  to  a 
Li  atom  and  a  H~  atom.  The  AGP  reduces  the  importance  of  this  unphysi- 
cal  configuration  by  making  the  second  and  third  largest  coefficients 
^/ery   small.  When  the  ratio  of  the  largest  coefficient  to  the  next 
largest  coefficient  approaches  infinity,  the  AGP  wavefunction  approaches 
a  generalized  AGP  wavefunction  (Weiner  and  Goscinski,  1981). 

The  AGP  at  the  equilibrium  internuclear  distance  provides  slightly 
less  than  half  of  the  correlation  energy  (see  Table  2-4).  This  result 
is  similar  to  the  results  obtained  by  other  investigators  (Jensen,  1980; 
Bessis  et  al.,  1967  and  1969;  Linderberg,  1980;  Kurtz  et  al . ,  1981). 
At  the  internuclear  distance  of  20.0  bohr,  however,  the  percentage  of  the 
correlation  energy  rises  to  92.  What  is  of  greatest  importance  is  the 
treatment  of  total  energy  differences  along  the  potential  curve  (see 
Table  2-5).  The  results  from  AGP  give  a  curve  that  is  similar  to  the 
full  CI  curve.  The  AGP  wavefunction  predicts  dissociation  into 
uncharged  atoms,  unlike  restricted  Hartree-Fock.  This  dissociation 
property  makes  the  AGP  a  useful  ground  state  for  the  investigation  of 
potential  energy  surfaces. 

2-8  The  random  phase  approximation  of  the  polarization  propagator 

The  AGP  is  now  used  as  a  ground  state  for  the  determination  of  exci- 
tation energies  with  the  random  phase  approximation  of  the  polarization 


67 


68 


o  o 


00    S- 


1/1 

1 ) 

<D 

eu 

U 

Q. 

c 

1/1 

CI) 

Cl) 

s. 

ec 

at\ 


69 


jropagator.  If  the  operator  manifold  is  transformed  so  that 


qijT  =  gVi+aj  -  gV/ai  ' 


(2-48; 


q  .  .  =  q  .a  .Ta .  -  q  .a  .  a  .  , 


(2-491 


where  j  =  j  with  opposite  spin,  then  the  polarization  propagator  will  be 
block  diagonal , 


p(e: 


.t._1\ 


0 


0 


«q;q». 


(2-50) 


for  a  self-consistent  RPA  polarization  propagator  (Weiner  and  Goscinski, 
1980a  and  1980b).  The  excitation  energies  are  obtained  from 


t  t    t      t 

«g  ;g  »r  gT  =  eQT 


(2-51) 


While  the  AGP  that  has  been  calculated  for  use  in  this  calculation  has 
been  variationally  optimized,  there  is  no  guarantee  that  it  is  the  AGP 
that  is  the  consistent  ground  state  of  the  RPA.  The  previously  cited 
proofs  only  state  that  the  consistent  ground  state  is  of  the  AGP  form. 
Nonetheless,  if  the  upper  block  of  equation  2-50  is  constructed  and 
diagonalized,  the  eigenvectors  form  excitation  operators  and  their 
adjoints  destroy  the  AGP  ground  state.  To  this  extent,  the  variationally 
determined  AGP  is  a  consistent  ground  state.  Of  course,  the  other  blocks 
of  the  propagator  have  been  neglected  and  are,  in  general,  not  equal  to 
zero. 
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2.9  Polarization  propagator  calculations  on  lithium  hydride 

The  lowest  excitation  energies  and  the  corresponding  excitation 
operators  are  calculated  according  to  the  previously  discussed  procedure 
(see  equation  2-51).  In  Table  2-6,  they  are  compared  to  full  CI  excita- 
tion energies.  Near  the  equilibrium  internuclear  distance,  the  agreement 
between  the  AGP-RPA  calculation  and  full  CI  is  good.  As  the  internuclear 
distance  is  increased,  the  results  grow  worse.  This  perhaps  is  due  to 
the  need  for  elements  of  f,   in  the  operator  manifold  since  more  than  one 
configuration  is  important  in  the  ground  state.  Specifically,  near  the 
equilibrium  internuclear  distance,  the  single  determinant,  molecular 
orbital  description  of  electronic  structure  is  qualitatively  correct. 
Excitations  are  approximately  single  replacements  from  a  reference  deter- 
minant. As  the  internuclear  distance  grows,  configuration  interaction 
in  the  ground  state  becomes  important.  The  two  configuration  wavefunc- 

tion  that  is  appropriate  for  such  regions  of  the  potential  energy  curve, 

2  2    2  2 
°l  aZ     "  al  a3  '  can  ^e  aPProx"imately  written  as  a  valence  bond  wavefunc- 

tion, 

det{(})Li)ls(l)a(l)0Lijls(2)3(2)4>Lij2s(3)a(3)4,H5ls(4)6(4)}  +    (2-52) 
det{<j>Lias(l)a(l)(j)Li)ls(2)3(2)(j>HJs(3)a(3)(j)Uj2s(4)3(4)} 

and  the  lowest  excitation  is  easily  characterized  as 

[\i,2s*H,ls  +  *H,ls*Li,2sJ  -"  *H,ls*H,ls  <2"53) 

which  is  a  single  excitation  from  the  valence  bond  wavefunction.  In  the 
delocalized  molecular  orbital  treatment,  the  same  excitation  must  be 
written  as 
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2  2     22     22     22     2        2 

Oj  02  "  CT^  °3  """ >  a2  a2  +  al  a3  "  al  a2a3  "  al  °3a2       (2-54) 


where 

ai  =  ♦Li.is  (2-55; 

a2     =     *L1,2s  +$H,ls  (2-56) 

a3     =     *Li,2s   "  *H,ls    .  (2"57) 

This  change  in  the  nature  of  the  electron  distribution  suggests  that 
perhaps  a  transformation  of  the  basis  of  excitation  operators  might  be 
useful  in  cases  where  more  than  one  configuration  is  dominant  in  the  CI 
wavef unction. 


APPENDIX 
CORRELATION  SELF-ENERGIES 

The  following  terms  should  be  excluded  from  the  self-energy  in  the 
calculation  of  electron  binding  energies  when  the  ASCF  plus  correlation 
self-energy  formulas  are  used.  These  exclusions  apply  only  to  the 
diagonal  elements  of  the  self-energies.  The  indices  refer  to  spin- 
orbitals.  Indices  a,b,c,...  refer  to  occupied  spin-orbitals,  and  indices 
p,q,r,...  refer  to  unoccupied  spin-orbitals. 

The  second  order  correlation  self-energy  is 


(a-i; 

Yc>2m    I  Y   Y  Y  <1'p1  [abxabj  |  ip>   1  r  r   r  <ia[  |pqxpq|  |ia> 

L         \c>  9      k         II  p+p    -c    -p  oil         I  F+p    -p    -p 

ii       c   a/i  b^i  p   t+£p  ea  %  L   a  p£i  q/i    t+ea  ep  £q 


The  diagonal  2p-h  TDA  correlation  self-energy  is 

(a-2; 


C.DIAG  2p-h  TDA  ,  .    ,  ,    .        .  ,  ,  . 

Y  (E)     =    I    y       y     y  <ip||ab><ab|hp> 

2  a£i   b£i   \  E+ep-£a-£b^<abHa^-<P 
1_  r     r       r  <iaj  |  pqxpq  |  |  ia> 


ii  2  a£.   ^.    I  E+ep-ea-cb+%<ab | | ab>-<pa | | pa>-<pb | | pb> 


l_  r      r         r  s  Id  M  pq^pq  II    d^ ^__ 

2  I  p^.   q£. E+ea-ep-eq-Js<pq 1 1 pq>+<aq 1 1 aq>+<ap 1 1 ap> 
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The  third  order  correlation  self-energy  is 


C,3  ,  <ia|  Ipqxpql  Irsxrs  I  I  ia>fl-6     6    .] 

ii  a  p  q  r  s  ^  £a  er  es^L+Vep     qJ 

,  <ip|  labxabl  Icdxcdl  |ip>fl-6     6. 

1  y  v  y  v  y         '  '  "        .      '  ac  1? 

'  4  £  £  t  h  k  (E+e„-e   -eJ[E+e  -e  -e. 


(A-3! 


abcdp  ^VVedJ  ^VVV 

<ib|  |rqxra|  |pbxpq|  |ia>f(l-6  .  6.    )(1-S     6-    )+(l-6    .)) 

Y  y  V  V  V  M  M        ^        ab  ia; v       pr  ip;    K       qi'J 

'Hill  i^a-p-qH^b-eq-er) 

„  <iq|  Icbxcpl  |aqxab|  |1p>((l-6ac6ia)(l-6pq6ip)  +  (l-6b.)) 

I  6  C  p  q  f^Wi^VT 

+  l  y  y  y  y  y  <ic[  labxabl  Ipgxpgl  |ic> 
+  I  v  v  y  V  Y  <ia|  Ipqxpql  [bcxbcj  [ia> 

4 1 :tlc  pi  K:vvK+vv£qJ 

+  I  v  V  V  V  V  5J  P 1  1  abxab  I  |  qrxqr  |  |  i  p> 

+  I  y  y  y  y  y  <ip|  |grxqr|  labxabl  |ip> 
4  a  11  q  r   ^p^a-SjUa^b-VSJ 

<ir|  |aq><ab|  |prxpq|  |  ib>fl-6   .] 

y  y  y  y  y       '  "_         "      ^     gv 

abpji       lE+ W£qJ  ^a+£b-£p-eri 

<ib|  |pr><pq|  |abxar|  |iq>fl-6   .1 

-  y  y  y  y  y  , ,  ri 

a  b  p  q  r  lE+ WerJ  ^W? 

<iq|  |acxab|  |pqxpc|  |ib>fl-5   .] 

Y  y  y  y  y  ci; 

abcpq         lE+VVeCJ  ^a+£b-£p-£qi 

ic|  |qbxqp|  |ac><ab|  |  i  p>  (l-6.-) 


inn-  ,,,,, .,  ^c ^  . 

abcpq         ^+ep  ea  %>  ^ea  ec  ep  eq 


E+en-£,-ehJl£,+er 
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1  V   V  V  Y   V    .<1'r  1  I  "i  p >< a b  |  |  rqxpq  |  |  ab> 


abpqr   lea+ W£qJ lea+eb"eq"eH 


1_  v  v  V  V   V    ,<"'al  1  ic><cb    1  pqxpq  |  ]ab> 

~  2  L  l  l  L  L    (e  +p'-e  -e       e. +e   -e   -e   i 

a  b  c  p  q    l  a  bb     p     qJ  l  b  ^c     p     qj 

.    1   r  v  v  v  v  <ip|  liaxabl  Iqrxqrl  |pb> 

2  i  i  I       U  +e  -e  -el (e  -e    1 — 

abpqr       ^ea  £b  £q  V  ^ea  V 

1  v  V  V  V  V  <ipl  I  ia><bc[  Ipqxaql  [bc> 

1  y  y  y  y  y  <j  a  |  |  j  r><rb  1  1  pqxpq  |  |  ab> 
+  2  L  L   L  L   I        le  +e. -e  -e    Me   -el 
abpqr       «■  a     b     p     q;  w     rJ 

1_  y  y  y  y  y  <icj  [ipxpq|  | ab><ab  1  [cq> 
"  2  a  b  C  p  q       ^VVVV^VV 
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